Short time dynamics with initial correlations 
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The short-time dynamics of correlated systems is strongly influenced by initial correlations giving 
rise to an additional collision integral in the non-Markovian kinetic equation. Exact cancellation of 
the two integrals is found if the initial state is thermal equilibrium which is an important consistency 
criterion. Analytical results are given for the time evolution of the correlation energy which are 
confirmed by comparisons with molecular dynamics simulations (MD). 



Although the Boltzmann kinetic equation is success- 
fully applied to many problems in transport theory, it 
has some serious shortcomings Among these, the 
Boltzmann equation cannot be used on short time scales, 
where memory effects are important |^,^ . In such situa- 
tions, a frequently used non-Markovian kinetic equation 
is the so-called Levinson equation Q,^. One remarkable 
feature of this equation is that it describes the formation 
of correlations in good agreement with molecular dynam- 
ics simulations Nevertheless, the Levinson equation is 
incomplete by two reasons: (i) It does not include corre- 
lated initial states, (ii) When the evolution of the system 
starts from the equilibrium state, the collision integral 
does not vanish, but gives rise to spurious time evolution. 
The latter point has been addressed by Lee et al. [Q who 
clearly show that from initial correlations there must ap- 
pear terms in the kinetic equation which ensure that the 
collision integral vanishes in thermal equilibrium. 

The aim of this letter is to derive the contributions 
from initial correlations to the non-Markovian Levinson 
equation within perturbation theory. We will restrict 
ourselves to the Born approximation which allows us to 
present the most straightforward derivation. The inclu- 
sion of higher order correlations can be found in |^,|| . 
The effect of initial correlations becomes particularly 
transparent from our analytical results which may also 
serve as a bench-mark for numerical simulations. 

The outline of this Letter is as follows. First we give 
the general scheme of inclusion of initial correlations into 
the Kadanoff and Baym equations in terms of the density 
fluctuation function. We show that initial correlations 
enter the kinetic equation as self-energy corrections and 
meanfield-like contributions in terms of the initial two- 
particle correlation function. An analytical expression 
for the time dependent correlation energy of a high tem- 
perature plasma is presented and then compared with 
molecular dynamics simulations. 

To describe density fluctuations we start with the 
causal density-density correlation function 



i(12l'2') = G2(121'2') - G(11')G'(22'), 



(1) 



where 1 denotes cumulative variables {xi,ti,..). 
G(l,2) = i(T*(l)*(2)+) and G2(121'2') are the one- 



and two-particle causal Green's functions. Their dynam- 
ics follows the Martin- Schwinger hierarchy 
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G2(121'2') = 

(5(1 - l')G(2, 2') - <5(1 - 2')G(2, 1') 

d3y(l,3)G3(l,2,3,l',2',3+), 



(2) 



where V{1,2) is the interaction amplitude and ^//(l) — 
J d2F(l,2)G(2,2+) is the Hartree self-energy. 
Using for G3 the polarization approximation 

G3(123l'2'3') = G(11')G(22')G(33') + 
G(11')L(232'3') + G(22')L(131'3') + G(33')i(12l'2'), (3) 

leads to a closed equation for L which is conveniently 
rewritten as integral equation 

L{1, 2, 1', 2') = Lo{l, 2, 1', 2') - Gh{1, 2')G{2, l') 

+ / d4Gif(l,4)G(4,l') / d31/(4,3)L(2,3,2',3+) (4) 



where (G§) ^ denotes the l.h.s. of the first equation (||) 
and we have taken into account the boundary condition 



(5) 



In the case that all times in (^) approach to, the right- 
hand side vanishes except Lq which represents, therefore, 
the contribution from initial correlations. They propa- 
gate in time according to the solution of (j^) Q 

Lo(121'2') = J dxidx2dx[dx2G^{l,xito)G'^{2,X2to) 

X Loo(xi , X2 , x[ ,x'2,to)Gi^(l', x[h)Gjj (2', x'^t^). (6) 

Here Lqo is the initial two-particle correlation function. 

Inserting (|^) into the first equation of (||) and restrict- 
ing to the Born approximation we obtain for the causal 
function [G^i^(l,2) = G^i(l,2) + y(l,2)G<(l,2)] 
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GffUl' 3)G(3, 2) = S{1 - 2) + 5i„it(l, 2) 

+ yd4{Eo(l,4) + I](l,4)}G(4,2), (7) 
c 

where the integration is performed along the Keldysh 
contour C with the self-energy in Born approximation 

E(l,2) = 

j d3d5y(l,3)Gff(l,2)F(2,5)Gff(3,5++)G(5,3+). (8) 

Two new terms appear due to initial correlations 

Eo(l,2)= J d3d5F(l,3)G//(l,2)T/(2,5)Lo(3,5,3+,5++), 

5init(l, 2) = y d3V{l, 3)Lo(l, 3, 2, 3+). (9) 

The integral form of (|^) is given in figure ^ from which 
the definitions M) are obvious. 



FIG. 1. The Dyson equation including density fluctuation 
up to second Born approximation. Besides the initial corre- 
lation term Sinit discussed in |^^, a new type of self energy 
So appears which is induced by initial correlations. Since the 
latter one contains interaction by itself, this term is of next 
order Born approximation. 

The equation for the retarded Green's function 
G^(l,2) = -ie(ii - t2)(G>{l,2) -f G<(1,2)), where 
G<(1,2) = {^+(2)^(1)) and G>(1,2) = (M'(l)«'+(2)), 
is derived from (0) as 

(G^V - - S^)G^ = ^(1 - 2) + S.!^,Al, 2) (10) 

and leads to the Kadanoff-Baym equation 

(G]^V-S^)G<-G<(G^V-S^) 

= (S + So)<G^ - G«(S + S]o)< + 5init - SU. (11) 

Using the generalized Kadanoff-Baym ansatz 

G<{h,t2) ^ lG''{h,t2)p{t2) - ip{tl)G^{h,t2) (12) 

we obtain the kinetic equation for the reduced density 
matrix p{t) = G< {t, t) 

d 

—p{k,t)=I{k,t)+jQ{k,t)+Ji{k,t) (13) 

with 



J(k,t) = -^Re / dti [ 
^ ' J J (2^ 



dqdp 
Kf 



V\q) 



xG"it,ti,k- q)G^{ti,t,k)G^{t,ti,p + q)G^{ti,t,p) 
X [pih, k - q)p{ti,p + q)il - p{ti,p)){l - p{ti, k)) 
-p{h,k)p{h,p){l - p{h,p + q)){l - p{h,k - q))] , (14) 
dqdp 



Ioik,t) = |lm J 



{2-Kny 



lV{q) 



xG'^it, h, k ~ q)G^it, to,k)G'^{t, h,p + q)G^{t, to,p) 
x{P-^ + q\Loo{p + k,to)\^), (15) 



I^{k,t) = Arc J dh J j^C^{q,tM)V\q) 

to 

xG'^it, h,k- q)G^{h,t, k) [p{h,k -q)- p{h,k)] (16) 



X1+X2 



where Loo(xi, X2, xa, X4) =< xi - a;2|Loo( 
^£^)|a:3 - a;4 > and Co{q,t,t') = / dxe""?^ < 
||Zyo(0)|| >. We hke to note that the equation (^ is 
valid up to second order gradient expansion in the spa- 
tial coordinate. This variable has to be added simply in 
all functions and on the left side of ( [l3| ) the standard 
mcanfield drift appears. 

The first part (|lj) is just the precursor of the Levin- 
son equation in second Born approximation ^ V^. The 
term ( [l^ ) coming from Eq leads to corrections to the 
third Born approximation since it is ~ y^£o- A more 
general discussion of higher-order correlation contribu- 
tion within the T-matrix approximation can be found in 
P,p^ and of general initial conditions in |jl^ . The second 
part (|l5|) following from S gives just the correction to the 
Levinson equation, which will guarantee the cancellation 
of the collision integral for an equilibrium initial state. 
Recently the analogous term in the collision integral has 
been derived by other means 

Multiplying the kinetic equation equation ( p"3[ - p^ ) with 
a momentum function 4>{k) and integrating over k, one 
derives the balance equations 



{m) = 



dk 



dk 



(27r?i)3'^' ' J {2TTnf 
For the standard collision integral follows 



0(fc)Xo. (17) 



{mi) 



xG^(t, h, k - q)G^(ti,t, fc)G^(i,ti,p -f g)G^(ti,t,p) 
xp{ti,k - q)p{ti,p + q){l - p(ti,p))(l - p{ti,k)) 



X <^ 0(fc) + ^{p) - ^{k -q)- <j>{p + q) 



(18) 



from which it is obvious that density and momentum 
{(f) — l,fc) are conserved, while a change of kinetic en- 
ergy (j) = k"^ /2m is induced which exactly compensates 
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the two-particle correlation energy and, therefore, assures 
total energy conservation of a correlated plasma . Ini- 
tial correlations, Eq. (p^), give rise to additional contri- 
butions to the balance equations MM- We get 



(0(fc)Xo) - ^ 



dkdqdp^^, . 




(27r?i)s 



q\Lo{p + k)] -^ „ ^ ) - c.c. 



2 " 2 

X {0(fc) + 0(p) - 0(fc - g) - 0(p + q)} (19) 

which keeps the density and momentum also unchanged 
and only a correlated energy is induced. The self-energy 
corrections from initial correlations which correct the 
next Born approximation, ([l^), would lead to 



rRe 



dkdq 
(27r?i)6 



dtiV\q)£o{q,t,ti)p{ti,k) 



xG^(t, ti,k^ q)G^{ti,t, A:)|0(fc - 9) - 



(20) 



which shows that the initial correlations induce a flux 
besides an energy in order to equilibrate the correlations 
imposed initially towards the correlations developed dur- 
ing dynamical evolution if higher than ^ correlations 
are considered. 

We will consider in the following only second Born ap- 
proximation ^ and have therefore to use from dTO) 



(21) 



G^(ii, i2, k) « -iQih - i2)e^^(*^-*i 
and for Lqq the first Born approximation 

— — \Loa{k+p)\ — 



= -7^^0(9) {po(fc)po(p)(l - Po(fc - q)){l - poip + q)) 
Ae 

-(1 - po(fc))(l - Po{p))iPo{k - q)){po{p + q))} . (22) 



where V denotes the principal value, Ae = ^ + 2^ — 



and po the initial Wigner distribution. 



Then the explicit collision integral (|T^ reads 



I{k,t) = — / dti 

fi Jto 



dqdp 



V\q) 



X COS 



2to 



P 

2™ 



(27r?i)6 

{k-qf ip + qr\{t-h) 



2m 2m / h 

X {piti,k- q)p{ti,p + q){l - p{ti,p) - p{ti,k)) 
X -p{h,k)p{ti,p){l - pih,p + q) - p{h,k - q))} (23) 

and the new term due to initial correlations dlq) is 
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dti 



dqdp 



V{q)V^[q) 



2m 



P 

2m 



{k-qf {p + q?\{t-h) 



2m 2m J fi 

X {pt){k - q)po{p + g)(l - Po{p) - Po{k)) 
X -po(fc)po(p)(l - Po{p + q)- Po{k - q))} ■ (24) 

To show the interplay between collisions and correla- 
tions, we have calculated the initial two-particle correla- 
tion function in the ensemble, where the dynamical inter- 
action V{q) is replaced by some arbitrary function VQ{q). 
Therefore the initial state deviates from thermal equilib- 
rium except when V{q) = Vo{q) and g(to) = go- 

The additional collision term, Xqj cancels exactly the 
Levinson collision term in the case that we have initially 
the same interaction as during the dynamical evolution 
(Vo = V) and if the system starts from the equilibrium 
p{t) = Po- Therefore we have completed our task and 
derived a correction of the Levinson equation which en- 
sures the cancellation of the collision integral in thermal 
equilibrium p5[. 




(Opt 

FIG. 2. The formation of correlation energy 
-fcorr = Etotai — finit " fcoii = fkin in & plasHia with De- 
bye interaction Vi. The upper panel compares the analytical 
results (^^ with MD simulations from for three different 
ratios of k_d to the inverse Debye length x = kd/k. In the 
lower panel we compare theoretical predictions for the inclu- 
sion of Debye initial correlations characterized by xq = acq/^: 
where x = kd/k ~ 1. 

On very short time scales we can neglect the change in 
the distribution function. Assuming a Maxwellian initial 
distribution with temperature T and neglecting degener- 
acy, we can calculate explicitly the collision integrals and 
obtain analytical results. We choose as a model interac- 
tion a Debye potential Vi{q) — 47re^?i^/[g^ -I- fi^Kf] with 
fixed parameter Ki = kjj and for the initial correlations 
Ki — kq. We obtain for the change of kinetic energy on 
short times from (^8|) and (p^ 
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§-^E^in{t) - S[Viqfm £[Voiq)V{q)]{t), (25) 
which can be integrated ||] to yield 

Ekin{t) = -Etotal — Einit{t) ~ Ecoll{t). (26) 

For the classical limit we obtain explicitly the time de- 
pendent kinetic energy 



nT 



Ax 



(27) 



where J-{y) = 1 — e** erfc(2;), r ~ tujp/\/2, x = kd/k 
and = Ane'^n/T — uj^T/m. The plasma parameter is 

given as usually by F = where = (jf;^)^^'^ is the 
Wigner-Seitz radius. 

In Fig. p, upper panel, we compare the analytical re- 
sults of (I?) with MD simulations using the Debye 
potential Vi as bare interaction. The evolution of ki- 
netic energy is shown for three different ratios x. The 
agreement between theory and simulations is quite satis- 
factory, in particular, the short time behavior for x = 2. 
The stronger initial increase of kinetic energy observed 
in the simulations at a; = 1 may be due to the finite size 
of the simulation box which could more and more affect 
the results for increasing range of the interaction. 

Now we include the initial correlations choosing the 
equilibrium expression (p2|) which leads to 



£^init if) 

nT 



2 ^iXq X 



[xT{xt) - XoT{xot)] , (28) 



where xq = kq/k characterizing the strength of the ini- 
tial Debye correlations (22) with the Debye potential Vb 
which containes kq instead of kd ■ Besides the kinetic en- 
ergy ( |2^ ) from initial correlations, the total energy -Btotai 
( p6| ) now includes the initial correlation energy which can 
be calculated from the long time limit of (E^ leading to 



£;total V3F3/2 



nT 



2{x + xq) 



(29) 



The result ( p6|) is seen in Fig. ||, lower panel. We ob- 
serve that if the initial correlation is characterized by a 
potential range larger than the Debye screening length, 
xq < 1, the initial state is over-correlated, and the cor- 
relation energy starts at a higher absolute value than 
without initial correlations relaxing towards the correct 
equilibrium value. If, instead, a;o = 1 no change of cor- 
relation energy is observed, as expected. Similar trends 
have been observed in numerical solutions ^ . 

In summary, in this Letter initial correlations are inves- 
tigated within kinetic theory. Explicit correction terms 
appear on every level of perturbation theory correcting 
the non-Markovian kinetic equation properly in a way 
that the collision integral vanishes if the evolution starts 
from a correlated equilibrium state. Furthermore, the 



conservation laws of a correlated plasma are proven in- 
cluding the contributions from initial correlations. It is 
shown that besides the appearance of correlation energy 
a correlated flux appears if higher than Born correlations 
are considered. 

Deriving analytical formulas for high temperature plas- 
mas allowed us to investigate the time dependent forma- 
tion of the correlation energy and the decay of initial 
correlations. The comparison with molecular dynamics 
simulations is found to be satisfactorily. Including ini- 
tial correlations the cases of over- and under-correlated 
initial states are discussed. While starting from equilib- 
rium the correlation energy does not change, for over- 
and under-correlated states the equilibrium value is ap- 
proached after a time of the order of the inverse plasma 
frequency. 
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